SRC  TR  87-156-rl 


On  Robust  Continuous-Time 
Discrimination 


By 

E.  Geraniotis 
and 

Y.  A.  Chau 


Report  Documentation  Page 

Form  Approved 

OMB  No.  0704-0188 

Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington 

VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  a  penalty  for  failing  to  comply  with  a  collection  of  information  if  it 
does  not  display  a  currently  valid  OMB  control  number. 

1.  REPORT  DATE 

1987 

2.  REPORT  TYPE 

3.  DATES  COVERED 

00-00-1987  to  00-00-1987 

4.  TITLE  AND  SUBTITLE 

5a.  CONTRACT  NUMBER 

On  RObust  Continuous-Time  Discrimination 

5b.  GRANT  NUMBER 

5c.  PROGRAM  ELEMENT  NUMBER 

6.  AUTHOR(S) 

5d.  PROIECT  NUMBER 

5e.  TASK  NUMBER 

5f.  WORK  UNIT  NUMBER 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

University  of  Maryland, Electrical  Engineering  Department, College 

Park, MD, 20742 

8.  PERFORMING  ORGANIZATION 

REPORT  NUMBER 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

10.  SPONSOR/MONITOR'S  ACRONYM(S) 

11.  SPONSOR/MONITOR'S  REPORT 
NUMBER(S) 

12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  unlimited 

13.  SUPPLEMENTARY  NOTES 

14.  ABSTRACT 

see  report 

15.  SUBIECT  TERMS 

16.  SECURITY  CLASSIFICATION  OF: 

17.  LIMITATION  OF 
ABSTRACT 

18.  NUMBER 
OF  PAGES 

44 

19a.  NAME  OF 
RESPONSIBLE  PERSON 

a.  REPORT 

unclassified 

b.  ABSTRACT 

unclassified 

c.  THIS  PAGE 

unclassified 

Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std  Z39-18 


ON  ROBUST  CONTINUOUS-TIME  DISCRIMINATION 
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University  of  Maryland 
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ABSTRACT 

Target  discrimination  problems  are  modeled  as  the  testing  of  binary  hypotheses  character¬ 
ized  by  continuous-time  observations  which  (i)  consist  of  distinct  signals  in  additive  white 
Gaussian  noise,  or  (ii)  are  the  output  of  stochastic  dynamical  systems  driven  by  white  Gaussian 
noise,  and  in  both  cases  have  partially  known  statistics.  In  particular,  the  signals  in  the  first 
model,  the  parameters  of  the  dynamical  systems  in  the  second  model,  and  the  autocorrelation 
functions  of  the  noise  in  both  models  belong  to  one  of  the  following  distinct  uncertainty 
classes:  classes  determined  by  2-altemating  capacities  and  classes  with  minimum  or  maximum 
elements.  Robust  discrimination  tests  with  a  fixed  observation  interval  and  sequential  tests  are 
derived  whose  likelihood  ratios  depend  on  the  least-favorable  pairs  of  parameters  in  the 
aforementioned  uncertainty  classes  and  are  shown  to  have  an  acceptable  level  of  performance 
despite  the  uncertainty.  For  tests  with  a  fixed  observation  interval  the  performance  measures 
considered  are  the  actual  error  probabilities  and  the  Chemoff  upper  bounds  on  them;  the  latter 
are  shown  to  preserve  their  desirable  asymptotic  properties  in  the  presence  of  the  uncertainties. 
For  sequential  tests  the  performance  measures  are  the  error  probabilities  and  the  average 
required  length  of  the  observation  interval  under  each  hypothesis. 


This  research  was  supported  in  part  by  the  Office  of  Naval  Research  under  contract  N00014-86-K-0013  and  in  part  by  the 
Systems  Research  Center  at  the  University  of  Maryland,  College  Park,  through  the  National  Science  Foundation’s  Engineering 
Research  Centers  Program:  NSF  CDR  8803012. 


I.  INTRODUCTION 


In  contrast  to  discrete-time  robust  detection  problems  (refer  to  the  tutorial  [1]),  continuous¬ 
time  robust  detection  problems-although  involved  in  several  situations  of  practical  interest- 
have  not  received  sufficient  attention  in  the  literature.  The  work  of  [2]  on  minimax 
continuous-time  discrimination  from  Poisson  observations  constitutes  an  exception  to  this. 
Motivated  by  practical  situations  encountered  in  target  discrimination,  we  investigate  in 
this  paper  detection  problems  from  continuous-time  observations  with  statistical  uncer¬ 
tainty  and  develop  detectors  (discriminators)  with  either  fixed-observation-interval  or 
sequential  processing  which  are  immune  (robust)  to  the  uncertainty  in  the  observations. 

Two  types  of  problems  of  testing  binary  hypotheses  are  considered  here:  the  first  type 
pertains  to  observations  that  are  signals  in  additive  white  Gaussian  noise,  namely 

Vi  Hi  :  Yt  =  s,(t )  +  Nt\  i  =  0,l 

where  0  <  t  <  T,  Si(t)  represents  the  signal  under  hypothesis  Hi,  and  Nt  is  a  Gaussian 
noise  process  with  covariance  E[NtN, ]  =  cr(t)cr(s)6(t  -  s). 

The  other  type  pertains  to  observations  that  constitute  the  path  of  the  state  in  a 
dynamical  system  driven  by  white  Gaussian  noise  described  by  the  following  stochastic 
equations  under  the  two  hypotheses 

P2  Hi  :  Tt  = 

dXi  ~  Si(t)Xtdt  -f  Nt;  i  =  0, 1 

for  0  <  t  <  T,  where  Nt  is  as  described  above,  Xq  is  a  normal  process  Xo  Af(0,l), 
and  Si(t)  is  the  parameter  of  the  dynamical  system.  Throughout  this  paper  we  place 
time  variables  in  the  subscript  for  the  random  processes  and  in  the  argument  for  the 
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deterministic  functions  of  time. 


Both  of  the  above  formulations  are  motivated  from  problems  of  practical  interest  in 
radar  target  discrimination  or  target  identification.  In  such  problems,  the  output  of  the 
signal  processor  must  indicate  which  of  several  targets  is  present  or  to  discriminate  between 
an  actual  target  and  a  decoy  or  chaff  cloud.  Usually,  this  discrimination  is  to  be  performed 
after  the  initial  decision  that  some  object  is  present  has  taken  place.  As  a  special  case 
of  the  two  models  above,  the  radar  detection  of  the  presence  of  a  target  could  also  be 
considered;  in  this  case  the  signal  under  the  null  hypothesis  is  identically  equal  to  zero. 
The  continuous- time  fixed-observation-interval  or  sequential  processing  considered  here 
can  model  the  processing  of  (i)  each  individual  pulse  of  the  radar  return  or  (ii)  of  several 
consecutive  pulses  of  the  radar  return. 

If  the  characteristics  of  the  channel  and  the  system  are  known,  then  Si(t)  and  o{i) 
(defined  above)  are  given.  In  reality,  however,  due  to  uncertainties  in  the  the  modeling  of 
the  targets  and  the  medium  through  which  the  radar  signal  and  the  radar  return  signal 
propagate,  the  waveforms  s;(t)  and  the  second-order  characteristics  of  the  additive  noise 
are  not  completely  known.  To  model  the  available  partial  knowledge  about  s,(f)  and  the 
noise,  we  consider  structured  uncertainty  classes  and  design  minimax  robust  discriminators 
for  problems  V\  and  TV 

We  consider  two  distinct  types  of  uncertainty  classes:  (i)  classes  determined  by  two  al¬ 
ternating  capacities  and  (ii)  classes  with  minimum  or  maximum  elements.  The  former  un¬ 
certainty  models  have  been  very  popular  among  the  robust  statisticians  and  engineers  (see 
[l]-[6])  because  they  include  several  useful  models  of  uncertainty  such  as  c-contaminated 
models,  total  variation  models,  band  models,  and  p-point  models  and  result  in  closed  form 
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expressions  for  the  least-favorable  elements  in  the  class  for  a  variety  of  performance  mea¬ 
sures.  The  latter  uncertainty  models  also  describe  useful  classes  of  uncertainty  and  result 
in  convenient  closed  form  expressions. 

The  design  philosophy  that  we  pursue  is  that  of  minimax  robustness,  which  has  received 
considerable  attention  for  in  the  past  fifteen  years  (see  the  tutorial  [1]).  According  to 
it,  the  worst-case  operational  conditions  are  identified  with  regard  to  some  performnace 
criterion  of  the  decision  design  and  the  optimal  such  design  for  these  conditions  is  derived. 
Subsequently  this  decision  design  is  employed  independently  of  actual  conditions  (which 
are  not  known,  except  for  the  fact  that  they  belong  to  some  structured  uncertainty  classes) 
and  it  is  shown  that  it  achieves  desirable  performance  despite  the  uncertainty.  In  this  paper 
the  design  criteria  for  problem  V\  are  the  probabilities  of  false  alarm  and  miss.  For  problem 
Vi,  since  there  is  no  closed  form  for  the  error  probabilities  under  the  two  hypotheses,  we 
derive  the  Chernoff  bounds  on  the  error  probabilities,  identify  their  connection  with  the 
discrimination  measures  known  as  the  divergences,  and  use  the  divergences  as  the  design 
criterion. 

This  paper  is  organized  as  follows.  In  Section  II,  we  first  review  the  definition  and 
known  properties  of  the  uncertainty  classes  characterized  by  2-alternating  capacities  and 
derive  some  new  useful  results  for  our  problems.  Then  we  describe  the  uncertainty  classes 
with  minimum  and  maximum  elements.  Finally,  we  cite  the  likelihood  ratio  functions 
for  problems  V\  and  Vi.  In  Section  III,  we  derive  the  minimax  robust  fixed-observation- 
interval  (block)  detector  for  V\  under  the  probability  of  error  criterion;  we  also  show  that 
the  error  probabilities  of  the  robust  test  converge  to  zero  exponentially  as  the  duration  of 
the  observation  interval  increases,  despite  the  uncertianty.  In  Section  IV,  we  derive  the 
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minimax  robust  block  test  for  Vi  when  the  performance  measures  are  the  Chernoff  bounds 
on  the  error  probabilities.  In  Sections  III  and  IV  both  the  uncertainty  classes  determined 
by  2-alternating  capacities  and  uncertainty  classes  with  minimum  and  maximum  elements 
are  considered.  In  Section  V,  we  derive  the  minimax  robust  sequential  detector  for  a 
homogeneous  Vi~ior  which  the  signals  s,(<)  (t  =  0, 1)  and  the  noise  variance  a(t)  are  time 
invariant;  the  uncertainty  classes  considered  here  are  of  the  minimum/maximum-element 
type.  Finally,  in  Section  VI  the  results  of  the  paper  are  summarized  and  conclusions  are 
drawn. 
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II.  PRELIMINARIES 


In  this  section  we  describe  the  two  uncertainty  class  models  that  we  use  in  this  paper 
and  the  associated  results  of  minmax  robustness.  We  also  derive  present  the  likelihood- 
ratios  for  problems  V\  and  Vi- 

II. A  Uncertainty  Classes  Determined  by  2- Alternating  Capacities 

First  we  provide  a  brief  but  self-contained  review  of  known  results  from  the  theory  of 
2-alternating  capacities  and  the  associated  theory  of  robustness  of  [6].  This  review  is  nec¬ 
essary  to  introduce  the  concepts,  notation,  and  fundamental  results  of  minimax  robustness 
within  these  uncertainty  classes,  which  we  use  repeatedly  in  this  paper.  Then  we  derive  a 
new  result,  Proposition  1,  which  extends  the  results  of  [6]  and  is  used  extensively  in  the 
proofs  of  our  results  on  robust  discrimination  described  in  Sections  III,  IV,  and  V. 
Definition:  A  positive  finite  set  function  tiona  sample  s  space  ft  and  associated  a-field  T 
is  called  a  2-alternating  capacity,  if  it  is  increasing,  continuous  from  below,  continuous 
from  above  on  a  closed  set,  and  satisfies  the  conditions  v(0)  =  0  and  v(AuB)  +  v(AC\B)  < 
v(A)  +  t>(R);  the  set  function  v  may  not  be  a  measure. 

Suppose  now  that  M  is  the  class  of  measures  on  (ft,.?7)  and  m  G  M  is  a  measure. 
Consider  the  uncertainty  class  which  is  determined  by  the  2-alternating  capacity  t;  as 
follows 

Mv  =  {m  €  M|m(A)  <  v(A),  VA  €  T,  m(ft)  =  t>(ft)}  (1) 

Thus  all  measures  in  the  above  uncertainty  class  have  the  same  upper  set  function  v.  When 
ft  is  compact,  several  popular  uncertainty  models  like  e-contaminated  neighborhoods  [3], 
total  variation  neighborhoods  [3],  band  classes  [4],  and  p-point  classes  (5]  are  special  cases 
of  this  model.  The  complete  description  of  one  of  these  classes,  the  f-contaminated  mixture 
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model  will  be  provided  later  in  this  section.  Fundamental  properties  of  the  uncertainty 
model  (1)  have  been  studied  by  Huber  and  Strassen  [6]  which-  enable  the  robustification 
of  the  average  risk  criterion  to  these  uncertainties.  We  state  the  relevant  properties  as 
Lemma  1. 

Lemma  Is  Suppose  i?o  and  v\  are  2-alternating  capacities  on  (Cl,  F),  (M0,  Mi)  are  the 
uncertainty  classes  determined  by  (noi^i)  as  in  (1),  and  the  two  classes  do  not  overlap. 
Then  there  exists  a  Lebesgue-measurable  function  irv  :  Cl  — >  [0,  oo]  such  that  the  average 
(Bayes)  risk  is  minimized 

9v0({7rv  >  0})  +  t»i({sr„  <  0})  <  6v0(A)  +  i>i(j4c)  (2) 

for  all  A  £  T  and  9  >  0;  6  can  be  interpreted  as  the  ratio  of  the  prior  probabilities  of  the 
two  hypotheses  H0  and  Hi,  whereas  {tt^  >  6}  can  be  interpreted  as  the  likelihood  ratio 
test  for  v0  versus  t>i.  Clearly,  (2)  together  with  (1)  imply  that 

9m0({ 7rv  >  0})  +  mi({7r„  <  0})  <  ^({^v  >  0})  +  ^({^v  <  0})  <  9v0(A)  +  t>i(Ac)(*) 

for  all  m0  £  Mo,  mi  £  Mj,  and  A  £  T\  this  inequality  establishes  the  minimax  robust¬ 
ness  of  the  test  based  on  7r„.  Furthermore,  there  exist  measures  (mo,  mi)  in  Mo  X  Mi 
such  that 

mQ({7ru  >  0})  =  v0({5Tv  >  #})  >  m0({7rv  >  0})  (3) 

mi({x„  <  9})  =  Vi({ttv  <  6})  >  <  9})  (4) 

The  quantity  7 r„,  which  is  termed  the  Huber-Strassen  derivative  of  the  classes  Mo  X 
Mi,  is  a  version  of  dmi/dmo  and  is  unique  a.e.  [mo];  it  plays  the  role  of  the  worst-case 
likelihood  ratio  for  the  two  uncertainty  classes  Mo  and  Mj.  The  dominance  properties 
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(3)  -  (4)  establish  the  existence  of  measures  in  the  classes  Mo  and  Mi  that  achieve  the 
upper  values  provided  by  v0  and  vj  for  sets  of  the  form  {irv  <  and  their  complements. 
The  measures  (mo,  mi)  are  termed  the  least-favorable  measures  over  Mo  x  Mj. 
Example:  Consider  the  e-contaminated  mixture  uncertainty  classes  of  probability  mea¬ 
sures  described  by  [2] 


Mj  =  {mj  €  M| mj  =  (1  -  Cj)m°  +  ejfhj,rhj(Sl)  =  m°(fl)  =  1},  j  =  0, 1. 

which  are  determined  by  the  known  nominal  measures  mg  and  m®  and  the  degrees  of 
uncertainty  Cq  and  fi  (0  <  e;  <  1  for  i  =  0, 1);  the  unknown  probability  measures  rhj 
(j  =  0,1)  are  allowed  to  take  any  arbitrary  values.  This  uncertainty  class  is  appropriate  for 
modeling  situations  in  which  the  probability  measures  (or  pdfs)  governing  the  observations 
are  convex  combinations  of  known  probability  measures  (pdfs)  and  arbitrary  probability 
measures  (pdfs).  Then  the  associated  2-alternating  capacities  are  for  j  =  0,1 


Vj{A)  = 


(!  —  e)m®(A)  -f  tj, 


0, 


A  =  0 


and  the  least-favorable  distributions  are 


dmo/d\ 


(1  -  eo)^mg/dA,  dm°/dmg  <  co 
drrii/dX ,  Co  <  droi/dmg, 


dmi/dX  = 


(1  -  e^dmi/dX,  c\  <  dm^/dm g 


ci(l  —  €i)dmg/dA,  dmj/dmg  <  Ci, 
where  A  is  the  Lebesgue  measure  and  O<C1<C0<00  are  constants  such  that  m i(fl) 

rho(fi)  =  1.  The  Huber-Strassen  derivative  rv  has  the  form 


7ru  =  dmi  fdrho  = 


1  -  ci 

1  -  e0 


min{co,max(ci,dmi/<img)}. 
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which  consists  of  a  censored  version  of  the  nominal  likelihood-ratio  dm\/dmQ. 

The  following  proposition  extends  the  dominance  properties  (3)  and  (4)  of  [6]  to  more 
general  functionals  than  the  error  probabilities  of  the  robust  test: 

Proposition  1:  Suppose  that  the  measures  (mo,  mi)  on  (fl,.?7)  belong  to  Mo  x  Mi 
characterized  by  (1)  and  that  x  is  a  real  variable. 

(i)  If  one  of  the  following  situations  holds: 

(a)  both  g(nv)  and  h(x )  are  nonnegative,  increasing  functions  of  irv  (the  Huber-Strassen 
derivative)  and  x,  respectively; 

(b)  g{nv)  is  a  nonnegative,  decreasing  function  of  7T„  and  h(x )  is  a  nonpositive  and  in¬ 
creasing  function  of  x; 

(c)  g(wv)  is  a  nonpositive,  increasing  function  of  nv  and  h(x)  is  a  nonnegative  and  decreas¬ 
ing  function  of  x; 

(d)  both  g( 7r„)  and  h(x)  are  nonpositive  and  decreasing  functions  of  and  x,  respectively, 
then 

/  g(*rj(x))h(x)mo(dx)  <  /  g(nv(x))h(x)m0(dx)  (5) 

Jn  J  n 

/  g(nv{x))h(x)mi(dx)  >  /  g(7rv(x))h(x)m1(dx)  (6) 

Jn  Jn 

(ii)  If  one  of  the  following  situations  holds: 

(a)  both  g(xv)  and  h(x)  are  nonnegative,  decreasing  functions  of  nv  and  x,  respectively; 

(b)  g{ 7TW)  is  a  nonnegative,  increasing  function  of  7rv  and  h(x )  is  a  nonpositive,  decreasing 
function  of  x; 

(c)  g( 7rv)  is  a  nonpositive,  decreasing  function  of  and  h(x)  is  a  nonnegative  and  increas¬ 
ing  function  of  x; 

(d)  both  g(i r„)  and  /i(x)  are  nonpositive  and  increasing  functions  of  irv  and  x,  respectively, 
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then 


J^g(rv(x))h(x)mo(dx)  >  J^g(Trv(x))h(x)rho(dx) 
j^g(xv{x))h(x)mi(dx)  <  J^g(irv(x))h(x)rhi(dx) 


(7) 

(8) 


where  mo  and  rh\  are  singled  out  by  Lemma  1. 

Proof:  We  only  prove  (6)  of  part  (i)  for  situation  (a),  since  (6)  for  other  situations,  (5) 
and  part  (ii)  can  be  proved  in  a  similar  way.  Applying  the  transformation  z  —  irv(x)  we 
have 

f  g(7Tv(x))h(x)m1(dx )  =  f  g{z)h{i:~l{z))mi(dz) 

Jo  1 

for  t  —  0, 1,  where  7T"1  is  the  inverse  image  of  ttv.  In  our  notation,  mi(dx)  =  mi({x  <  X  < 
x  +  dx})  =  dM\(x),  where  the  associated  distribution  function  is  M\(x)  =  mi({X  <  x}). 
Now,  assuming  that  the  conditions  in  (i)  hold  (i.e.,  g'(x )  >  0,  h'(x)  >  0,  g(x )  >  0  and 
h(x)  >  0,  with  the  prime  denoting  a  derivative  with  respect  to  the  argument)  and  using 
the  above  transformation  and  integration  by  parts,  we  have 


f  g(-Kv(x))h(x)mi(dx)  -  \  g{irv(x))h(x)rhi(dx) 

J  o  Jo 

=  I  9{z)h{^~l{z))mi{dz)-  f  g(z)h(n~l(z))rhi(dz) 

J  n  Jo 

=  f  mi({Z  <  z})d[g{z)h(x~\z))]  -  [  mi({Z  <  z})d[g(z)h(n~\z))] 

Jo  Jo 

=  /  ( mx{{Z  <  z})  -  mi ({Z  <  z}))h(n~1(z))dg(z) 

Jo 

+  f  (mi({Z  <  z})  -  mi({Z  <  z}))g(z)dh(ir-\z)) 

J  o 

=  [  (rhi({7rv(X)  <  z})  -  mxCI^X)  <  z}))h(nZ \z))g\z)dz 

Jo 

+  f  (mrCKPO  <  *„(<)})  -  mx({*„(*)  <  *v(t)}M*v(t))h'(t)dt  >  0 

Jo 

where  Z  —  7r„(X);  the  last  inequality  derives  from  (4)  since  z  =  x„(t)  >  0  for  all  t. 
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Remark  1:  If  either  g(x)  =  1  or  h(x)  =  1  for  all  x,  i.e.,  one  of  the  functions  g  or  x 
is  absent  from  the  integrands  of  (5)-(8),  the  inequalities  in  (5)-(8)  still  hold;  in  this  case, 
the  nonnegativity  or  nonpositivity  of  the  function  involved  is  not  a  necessary  condition. 

Remark  2:  Lemma  1  and  Proposition  1  hold  even  if  the  2-alternating  capacity  v0  is 
itself  a  measure.  In  this  case,  the  uncertainty  class  Mo  has  a  single  element  vq. 

II.B  Uncertainty  Classes  with  Minimum  or  Maximum  Elements 
In  our  cases  of  interest  these  uncertainty  classes  are  of  the  form 

s0,u  =  {*o(0M0  <  o  <  t  <  r},  (9) 

Si.L  =  {*i(0M0  >  *il(0,  0  <t<T},  (10) 

and 

Eu  =  M0K0  <  oi/(0.  0  <  t  <  T}.  (11) 

Uncertainty  classes  (9)  and  (11)  are  characterized  by  maximum  elements  $u(t)  and  cru(t), 
respectively,  while  (10)  is  characterized  by  the  minimum  element  si(t).  These  uncertainty 
classes  are  easy  to  characterize  and  for  the  performance  measures  of  interest  their  least- 
favorable  elements  turn  out  to  be  the  maximum  and  minimum  elements  involved  in 
the  definition  of  the  classes.  These  classes  are  appropriate  in  situations  in  which  upper 
and/or  lower  bounds  on  the  range  of  the  system  and  noise  parameters  are  available  (can 
be  estimated).  To  guarantee  that  the  classes  S0,U  and  Sj  l  d°  n°t  overlap  we  impose 
the  condition: 

*il(0  >  *01/(0.  o  <  t  <  T  (12) 
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II. C  Likelihood  Ratios  for  Problems  V\  and  Vi 

In  the  following  we  derive  the  likelihood  ratio  functions  for  V\  and  Vi-  Notice  that 
the  observation  process  in  V\  can  be  written  in  differential  form  as  follows: 


Hi  :  dYt  =  Si(t)dt  +  dWt  ,  t  =  0, 1;  0  <  t  <  T 


where  Wt  is  a  Wiener  process  with  zero  mean  function  and  intensity  function  a(t).  Using 
this  differential  form  of  observations  and  manipulations  analogous  to  the  derivation  of 
generalized  likelihood  ratio  (see  [7]  and  [8]),  we  give  Lemma  2  without  proof  for  the 
likelihood  ratio  function  of  Vi- 

Lemma  2:  The  likelihood  ratio  function  for  V\  has  the  following  form 


dPo  Jo  K0] 


2K0]2 


(13) 


where  Pi  for  i  =  0, 1  are  the  joint  probabilities  of  the  observations  {Yt,  0  <  t  <  T}  under 

Hi- 

By  using  the  explicit  form  of  the  solution  for  Xt  (and  thus  for  Yt )  to  the  stochastic 
differential  equation  in  Vi  (see  Section  4.4  of  [9]),  we  can  derive  the  likelihood  ratio  function 
of  Vi,  which  is  stated  here  as  another  lemma  and  its  proof  is  omitted. 

Lemma  3:  The  likelihood  ratio  function  for  Vi  has  the  following  form 


in  Lr  =  i„  §•  =  r  YldYl  _  [T  [»;«y  :  jfoW) \Y}dt 

dPo  Jo 


WOP 


1: 


2  WOP 


(14) 


where  P,  denotes  the  joint  probabilities  of  the  observations  {Yt,  0  <t<T}  under  Hi,  for 
i  =  0,l. 

The  likelihood  ratio  test  is  used  as  the  decision  rule,  since  it  minimizes  the  expected 
risk  and  provides  one  side  of  the  inequality  for  the  saddle-point  description  of  the  minimax 
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robust  scheme  [refer  to  equation  (*)].  The  criterion  used  for  deriving  the  minimax  robust 
scheme  for  V\  is  the  probability  of  false  alarm,  a(mo,  Lt),  and  the  probability  of  miss, 
P{tu\,Lt)-  Under  mismatch,  that  is  when  the  statistics  mo  and  mx  for  which  Lt  = 
dxh\ / TOq  is  derived  are  different  than  the  actual  too  and  TOj  governing  the  observations, 
the  probabilities  of  false  alarm  and  miss  become  a(mo,  Lt)  and  ^(ttii,Lt)-  For  V2  we  use 
the  Chernoff  bounds  on  a  and  (i  as  the  criterion  for  robustness. 

To  prove  minimax  robustness  one  has  to  show  that 

a(m0,  Lt)  <  <*(m0,  LT )  <  a(m0,  Gt)  (15) 

and 

/3(mi,Ir)  <  P(muLT)  <  P(rni,GT)  (16) 

over  all  to0  €  Mo,  toj  €  Mi,  and  arbitrary  decision  tests  Gt  based  on  the  observation 
interval  [0,Tj;  the  right-hand  inequalities  in  the  equations  above  are  satisfied  immediately 
once  the  likelihood-ratio  Lt  has  been  derived.  Therefore,  in  the  following  sections  we  only 
need  to  prove  the  left-hand  inequalities  in  (15)-(16)  for  the  various  cases  of  interest  in 
order  to  establish  the  minmax  robustness  of  the  likelihood  ratio  tests. 

For  convex  classes  Mo  and  Mi  the  above  inequalities  are  equivalent  to 

maxmoa(m0, Lt),  maxmi/3(mu LT)  (17) 

As  shown  in  [6]  the  least-favorable  pair  (mo,  mi)  for  general  classes  of  probability 
measures  Mo  and  Mi  of  the  form  of  (1)  for  the  problem 

maxa(mo,I),  max/3(m1,  L) 

tn  0  mi 

where  L  =  dmi/dmo,  is  also  the  least-favorable  pair  for  each  one  of  the  following  problems: 
(i)  rnax{maxmo  q(to0,  I),maxmi  /3(toi,I)}; 
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(ii)  maxmi  (3(m\,L)  subject  to  max,,,,,  a(m0,  L )  <  ao; 

(iii)  maxm0)mi(Aa(m0,  L)  4-  (1  -  \)P(mi,L))  with  A  the  known  prior  probability  of  H0. 
In  the  above  discussion  we  omitted  So,  Si  and  a  in  the  argument  of  a  and  (i.  However, 

when  we  consider  the  uncertainty  classes  (9)-(ll),  s0,  Si  and  a  are  the  variables  to  be 
robustified  and  the  above  statements  are  still  valid.  In  the  following,  we  will  show  explicitly 
in  our  notation  only  the  arguments  of  a  and  /?  that  are  to  be  robustified. 
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III.  MINIMAX  ROBUST  BLOCK  DISCRIMINATION  FOR  Vx 
III. A  Uncertainty  in  the  Signal 

1.  Uncertainty  Class  Characterized  by  2-AIternating  Capacities 

In  order  to  apply  Lemma  1  and  Proposition  1,  the  space  (II ,T,  m,)  t  =  0, 1  is  identified 
as  the  observation  interval  [0,T]  and  the  tr-field  T  is  the  one  generated  by  all  subsets  of 
[0,T].  In  this  subsection,  we  consider  the  particular  case  where  So(t)  is  a  known  nonnega¬ 
tive  time- varying  function;  thus  there  is  now  uncertainty  under  hypothesis  IIq.  Then  the 
class  So  contains  a  single  measure  and  is  defined  as 

S0  =  -Mmo(A)  =  s0(t)\(dt),  VA  €  T,  }  (18) 

where  A  G  F  and  A  is  the  Lebesgue  measure.  Therefore,  in  this  case,  the  probability  of 
false  alarm  remains  unchanged  over  Sq:  a(so(t),Lx)  =  a(so(t),Lx),  so  that  we  only  need 
to  consider  the  probability  of  miss.  The  uncertainty  class  under  hypothesis  Hi  is  defined 
as 

Sx  =  {sx(t)|mi(A)  =  f  |si(t)|A(df)  <  ri(A),  VA  G  T,  [  =  »i(ft)  =  P,}  (19) 

Ja  J  n 

This  definition  allows  for  Si(t)  to  take  both  positive  and  negative  values,  since  it  is  |si(t)! 
which  is  involved  in  the  definition  of  the  measure  mi  and  not  Si(t)  itself.  It  is  assumed 
that  Si(t)  /  so(t)  for  some  t  G  [0,T]  for  all  elements  of  the  class  (19)  so  that  there  is  no 
overlap  between  Si  and  So- 

As  mentioned  in  Section  II,  the  uncertainty  class  (19)  includes  several  useful  uncer¬ 
tainty  models  (including  four  of  the  most  popular  uncertainty  class  models).  In  adition  to 
those  classes,  the  following  two  other  uncertainty  classes  based  on  norms  can  be  considered: 

S2  =  {<si(0l  /  lsi(0  -  Si(t)\dt  <  ^,  VA  G  F)  ( l}norm )  (20) 

J  A 
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and 


S3  =  {51(f)!  sup  |$j(f)  —  si(OI  ^  ^  )  ( L°°norm )  (21) 

t 

We  can  show  that  Sj  C  Si  for  j  =  2,3  as  follows.  For  any  «i(t)  £  S2  we  have 

JA  |«i(0ld<  -  JA  l«i(0ld*  <  JA  l«i(0  -  S°(t)\dt  <  S, 

therefore,  if  we  define  v2(A)  =  fA  |s°(t)|d/  +  6,  then  s\ (t)  £  Si  as  well.  Similarly,  for  any 
Si(t)  €  S3  we  have 

f  |si(t)|dt  -  f  |sj(t)|dt  <  f  |«i(i)  —  <  sup  —  ^°(f)|A(J4)  <  6A(v4), 

Ja  Ja  Ja  t 

therefore,  if  we  define  vz(A)  =  fA  |si(f)|d<  +  £A(A),  where  A(A)  <  A  (Cl)  =  T,  then  Si(t)  £ 

Si,  as  well.  Subsequently,  since  Sj  C  Si  for  j  =  2,3,  the  following  inequality  holds 


max  /3(si(t),Xr)  <  max  (3(si(t),LT),  j  =  2,3. 

*i(t)eSj  «i(«)eSi 


(22) 


where  Ji(t)  =  |si(<)|  and  sQ(t)  are  the  least-favorable  densities  for  Si  vesus  So  and  Lt 
is  the  likelihhod  ratio  for  signals  Jj(t)  and  so(t)  obtained  from  (13).  Since  there  is  no 
uncertainty  under  Hq,  the  probability  of  false  alarm  trivially  satisfies 


max  o(«o(l),ir)  <  max  a(s0(t),LT),  j  -  2,3, 

*i(<)GSj  *i(<)Gsl 

with  equality  for  so(t)  =  So(<)- 

Proposition  2  provides  the  minimax  robust  test  for  signal  uncertainty  in  V\. 
Proposition  2:  For  problem  V\  and  signal  uncertainty  within  classes  (19)  and  (18),  if 
0  <  •si(f)  <  s0(t),  s'0(t)  <  0,  and  cr'(t)  >  0  for  0  <  t  <  T,  then  the  likelihood  ratio  test 
:  based  on 


In  Lt  =  In 


dPj 

dP0 


si(t)  -  s0(t) 

WOP 


Ytdt- 


fT  MW-M)? 

Jo  2[(j(t)]2 


(23) 
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where  Si(t)  and  «o(0  =  «o(0  are  the  least-favorable  elements  for  Si  vs  Sq  of  (19)  and  (18) 
singled  out  by  Lemma  1,  is  minimax  robust;  that  is,  the  two  types  of  error  probabilities 
for  this  test  satisfy 

a(s0(t),  Lt)  <  <*(so(0»£r) 

and 

P(si(t),LT)<l3(h(t),LT)  (24) 


Proof:  From  Lemma  2  we  can  obtain  (23)  under  worst-case  operational  conditions.  Since 
So(t)  is  known,  the  inequality  about  the  probability  of  false  alarm  is  trivially  satified  with 
equality.  Next  we  prove  (24).  Let 


ri 


-i 


T  [ii(0)2  -  MOP 
2W0P 


*  --i  n 


TfSl{t)  S°M}2[a(t)]2dt, 


K0]a 


and  define 


f  _  fTh(t)-s o(«)_/4\j,.  ;  _  rTh(t)-s0(t) 

h  ~  Jo  K<)]2  l(  )d  ’  h  ~  Jo  KOI2  li  )d 


From  ii(t)  <  so(t )  we  have 


h  “  Jo  kop  l( } 


KOI 

T  ^i(t)  -  s0(t) 


>  /  Mil 
-  Jo  [a 


K*)]2 


jsi(t)|dt. 


If  so(<)  7^  0,  then  since  a'(t)  >  0  and  s'0(t)  <  0  from  (c)  of  part  (i)  in  Proposition  1,  we 
use  (6)  to  obtain 


h  -  L  )|5,(<)l<“ 
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Finally,  since  f\  >  fi,  it  is  easily  observed  that 


^(t),Lr)  =  SI 

=  fKh(t),LT) 

since  $[(ln  tj  +  ri  —  «)/v/r2]  is  a  decreasing  function  of  x. 

2.  Uncertainty  Classes  with  Minimal  and  Maximal  Elements 

The  minimax  robust  discriminator  for  the  uncertainty  classes  (9)  and  (10)  is  given  by 
the  following  proposition. 

Proposition  3:  For  V\  and  the  uncertainty  classes  (9)  and  (10),  if  siz,(f)  >  sou(t)  for 
0  <  t  <  T,  then  the  likelihood  ratio  test  based  on 


<  $ 


In  r)  +  ri  -  fx 


v/^2 


m  ir  =  *  -  [T  MO zjMyldt  _  r  MSLzMZ*  (25) 
dPo  Jo  [cr(<)]2  Jo  2  [a{t)Y 

where  s0{t)  =  s0u(t )  and  sj(f)  =  sn(t),  is  minimax  robust,  that  is,  the  error  probabilities 
satisfy 

a(s0(t),LT)  <  a(s0(t),LT)  (26) 


and 

/3(*i  (0>£t)</W*)£t),  (27) 


Proof:  We  only  give  the  proof  of  (26),  since  (27)  can  be  proved  in  an  identical  manner. 
Let 


ri 


moi2  -  im? 

2K0P 


dt] 


r 2  = 


fT 

Jo  x  WOP 


}2w<)]2* 


and  define 


a-/ 


T  h(t)  -  s0(t) 

WOP 


-  rT  s 
s0(t)dt ;  fo=  ~ 
Jo 


Sl(t) 


WOl 
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i 


Then  fo  <  fo  and 

«(*>(<).  It)  =  p»(b^  >>«  ■»)  =  <? 

<  <3  — —  =«(«o(<)i^t) 

V^2  J 

where  we  used  the  fact  that  Q[(Irt]-\-  r\  —  x)fy/r^  is  an  increasing  function  of  x  and  that 
«i(0  =  slL(t)  >  Sou(t)  =  s0(t). 

Notice  that  the  distance  of  the  least-favorable  pair  |ii(t)  -  So(*)|  =  <si l(<)  -  So£/(t)  is 
the  shorest  of  all  distances  of  the  pair  (so(t),  si (t))  €  So^  x  Si  L-  This  result  is  intuitively 
pleasing,  since  the  closer  the  two  signals  (so(t),  s\(t))  are,  the  higher  the  error  probabilities 
are  expected  to  be. 

III.B  Uncertainty  in  the  Noise  Variance 

1.  Uncertainty  Classes  Characterized  by  2-Alternating  Capacities 
In  this  subsection,  we  consider  the  following  uncertainty  class  for  a(t) 

E  =  {<r(t)|m(A)  =  J  [<7(t)]2A(dt)  <  v(A),V/l  €  T,  J  [<r(t)]2A(cft)  =  v(ft)  =  p#}  (28) 

where  m(A)  defines  a  measure  on  (ft,  T)  with  ft  and  T  interpreted  as  in  III. A.  In  this  case 
Mj  =  X  and  Mo  =  {A}  and  we  assume  that  <r(t)  ^  1  for  some  t  £  [0,T]  such  that  there 
is  no  overlap  between  Mi  and  Mo-  Note  that  the  constraint  J^[<r(t)]2A(<fi)  =  v(ft)  =  ps 
implies  a  fixed  noise  power.  The  minimax  robust  discriminator  for  V\  is  characterized  by 
the  following  proposition. 

Proposition  4:  For  problem  V\  and  noise  uncertainty  within  the  class  (28),  if  (i)  s\(t)  > 
$o(t)  and  s[(t)  <  s'0(t)  or  (ii)  si(t)  <  sq(0  and  Sj(t)  >  s'0(t)  for  0  <  t  <  T,  then  the 


likelihood  ratio  test  based  on 


In  Lt  —  In  — -r-  =  / 

dP0  Jo 


dP\  fT  si{t)  -  s0{t) 


Ytdt-  ( 
Jo 


T  M01!  -  MOP 
2(*(0P 
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where  [ct(<)]2  *s  the  Huber-Strassen  derivative  for  Mi  =  S  vs  Mo  —  {A}  singled  out  by 
Lemma  1,  is  minimax  robust;  that  is,  the  error  probabilities  satisfy 


a(a(t),LT )  <  a(ar(t),LT)  (30) 

and 

P(<r(t),LT)  <  (3{cr{t),LT)  (31) 


for  all  thresholds  r/  satisfying 


Proof:  We  only  prove  the  part  with  assumption  (i),  since  under  (ii)  it  can  be  proved  in 
a  similar  way.  Let 


f  fT  5l(0  “  5°(0  /fUf  f  _  fT  *1(0  ^°(^).  o-  ( t\r1t 

/o  =  L  [»t01~^( 0  ’  h  "  Jo  HW  l{  >  ’ 


Under  the  assumption  (i)  [si(t)  -  so(0]2  is  a  decreasing  function  of  t,  thus  from  (a)  of  part 
(ii)  in  Proposition  1  we  use  (8)  to  obtain  r 2  <  tV  Then  (30)  follows  by  noticing  that 


a(a(t),LT)  =  P0(ln^>\nTi 

V  dp0  ) 

In  t/  +  n  -  /o 


=  Q 


In  t)  +  ri  -  Jo 


<  Q 


\/h 


y/ri 
a(cr(t),LT ) 


where  we  used  the  fact  that,  under  (32),  In  T)  +  ri  —  fo  >  0  and  thus  Q[(ln  rj  +  rj  fo)/x] 
is  an  increasing  function  of  x.  The  inequality  in  (31)  can  be  proved  in  a  similar  way  by 
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using  the  fact  that  r2  <  f2  and  that,  under  (32),  In  rj  +  rx  -  f\  <0,  which  in  turn  implies 
that  $[(ln»7  +  ri  -  fi)/x]  is  an  increasing  function  of  x. 

2.  Uncertainty  Classes  with  a  Maximal  Element 

The  minimax  robust  scheme  for  the  uncertainty  class  (11)  is  given  as  the  following 

proposition. 

Proposition  5:  For  problem  V\  and  the  uncertainty  class  (11)  the  likelihood  ratio  test 
based  on 


dP, 

In  Lt  =  In  — ir- 

dP0 


[ T  fl Q — ?^Ytdt 

Jo  KOI2 


h(t)l2-[so(t)]2 

2  m? 


dt 


where  a(t)  =  <?i /(<),  is  a  minimax  robust  test;  that  is 


(33) 


a(tr(f),l-r)  <  a(d(t),i/T) 


(34) 


and 


for  all  thresholds  r?  satisfying 


,,  , .  r  MO-MO]2,* 

2 w)f-dl 


(35) 


(36) 


Proof:  We  prove  (34)  and  (35)  following  steps  similar  to  those  in  the  proof  of  Proposition 
4  and  using  the  inequality 

fTrM0  -  MO 


T2 


-  (irlW4=f' 


III.C  Uncertainty  in  the  Signal  and  Noise  Variance 

In  this  subsection,  we  treat  the  general  case  where  the  signals  and  the  variance  of  the 
only  known  to  belong  to  the  uncertainty  classes  described  in  II.  By  examining 


noise  are 
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Propositions  2-5  we  notice  that  the  minimax  robust  scheme  for  this  general  case  can  be 
derived  directly  from  the  combination  of  their  results.  Let  us  now  address  the  minimax 
robust  tests  for  different  combinations  (different  uncertainty  classes  corresponding  to  the 
signal  and  the  noise)  as  four  corollaries.  Their  proofs  are  omitted,  since  they  can  be 
deduced  from  the  proofs  of  Propositions  2-5. 

Corollary  1:  For  problem  Vi  with  signal  and  noise  uncertainty  within  classes  (19)  and 
(28),  suppose  that  s0(t )  is  a  known  time- varying  function,  0  <  si(t)  <  so(t),  s\(t)  >  -So(t), 
s'0(t)  -  0,  and  cr'(t)  >  0  for  0  <  t  <  T,  where  sj(t)  and  cr(t)  are  the  least  favorable 
elements  of  classes  (19)  and  (28)  singled  out  by  Lemma  1  [applied  to  (19)  . ys  (18)  and  (28) 
vs  {A}],  then  the  likelihood  ratio  test  based  on 

dPi 


In  Lt  =  In  =  [ 
dP0  Jo 


T  s 


M‘)  -  *Myldt  _  fT  giCOf  zMjWdt 

Jo 


i<w  j„  mo? 

is  minimax  robust;  this  implies  that  the  error  probabilities  satisfy 

o(a(t),ir)  <  a(a(t),Lr) 


(37) 


and 


f)(si(t),cr(t),LT)  <  P(si{t),cr(t),LT ) 


(38) 


(39) 


for  all  thresholds  tj  satisfying 


llnj?l  <  / 

J  Q 


fsi(t)  -  So]2 


dt. 


(40) 


2[d(<)]2 

Corollary  2:  For  problem  V\  with  signal  and  noise  uncertainty  within  classes  (19)  and 
(11),  suppose  that  so(t)  is  known,  0  <  si(<)  <  so(t),  s'0(t )  <  0  o'(i)  >  0  for  0  <  t  <  T, 
where  s\(t)  is  the  least-favorable  element  of  (19)  [for  testing  (19)  vs  (18)]  singled  out  by 
Lemma  1  and  a(t)  —  au(t)>  then  the  likelihood  ratio  test  based  on  (37)  is  a  minimax  robust 
test;  that  is,  the  error  probabilities  satisfy  (38)  and  (39)  for  all  thresholds  t]  satisfying  (40) 
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Corollary  3:  For  problem  V\  with  signal  and  noise  uncertainty  within  classes  (9),  (10), 
and  (28),  if  su(t)  >  sou(i)  and  s'1L(t)  <  s^;(t)  for  0  <  t  <  T,  where  Jo(t)  =  s0i/(t), 
3j(t)  =  suit)  and  d(t)  is  singled  out  by  Lemma  1  [for  testing  (28)  vs  {A}),  then  the 
likelihood  ratio  test  based  on 


.» iT  =  =  fr  mYtdt  -  r  MtWf4  (41) 

dPo  h  [a(()]!  Jo  2[^(<)1 


is  minimax  robust;  in  this  case  the  error  probabilities  satisfy 


a(3o(t),o(t),lT)  <  a(ia{t),c(t),Lr) 


(42) 


and 

P(si(t),a(t),lT)  <  P(si(t),d(t),LT )  (43) 


for  all  thresholds  rj  satisfying 


[h{t)~s0{t)}2 

2  m? 


(44) 


Corollary  4:  For  problem  Vi  with  signal  and  noise  uncertainty  within  classes  (9)-(ll), 
suppose  s\i{t)  >  s0u(t)  for  0  <  t  <  T,  then  the  likelihood  ratio  test  based  on  (41)  with 
io(0  =  Soi/(<),5i(t)  =  sn(t)  and  &(t)  =  a {/(*)>  is  minimax  robust;  that  is,  (42)  and  (43) 
hold  for  all  thresholds  satisfying  (44). 
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III.D  Asymptotic  Performance 

It  is  interesting  to  examine  for  the  general  case  of  1II.C  the  behavior  of  the  two  types  of 
error  probability  as  T  (the  length  of  the  observation  interval)  increases.  It  is  particularly 
desirable  to  ascertain  that  the  robust  likelihood  ratio  tests  described  in  Corollary  1-4 
maintain  their  robustness  asymptotically  (for  large  T ).  Under  mismatch  and  for  the 
threshold  17,  the  Chernoff  bounds  on  the  two  error  probabilities  have  the  forms 

<*l/(-s,  s0(t ),  a(t),  Lt)  =  exp{-S77  -  TC0{s,  Lt)}  (45) 


and 


A/(s,  ~  exP(ST?  ~  TCi(s,  Lt)} 


(46) 


for  all  s  6  [0, 1].  In  (45)  and  (46),  the  Chernoff  distances  Cj(s,  Lt)  ( j  =  0, 1)  under 
mismatch  are  defined  as 

Co(s,Lt)  =  -HEo{Lt}]/T  (47) 


and 


C\U,l7)  =  -  In  {E1{lT))IT. 


(«) 


For  V\  we  can  easily  obtain  that  Ej{L‘}  are  given  by 

Ei{f}  =  *{«*  f  expHjf  («) 

for  j  =  0,1.  If  we  define  Zj-  =  Jg  then  Zt  is  normally  distributed  with 


mean 


l 


TSi(t)-s0(t)  m  . 

mw'  A) 


under  Hj,  and  variance 


23 


Upon  substitution  into  (47)-(49)  we  obtain 


C0(s,  lT)  =  ^[ci s2  +  c2s  -  s  f 

1  Jo 


Tii(t)-s0(t)  . 

li?  ~So{t)dt] 


and 


where 


Ci(s,Lt)  =  j;[cis*-c1s  +  sJv  (!)<<'] 

T  MOP  -  MOP 


dt. 


(50) 


(51) 


(52) 


[*(0P  Jlvvv,,",v*'y0  2[dr(t)]2 

The  following  proposition  guarantees  the  robustness  of  the  likelihood  ratio  tests  in 
Corollaries  1-4  for  the  Chernoff  bounds  on  the  error  probabilities. 

Proposition  6:  For  problem  V\  with  the  uncertainty  models  for  the  signal  and  the  noise 
described  in  Corollaries  1-4,  the  likelihood  ratio  tests  and  the  least  favorable  pairs  provided 
there  are  minimax  robust  and  least-favorable,  respectively,  for  the  Chernoff  bounds  on 
error  probabilities;  that  is,  the  following  inequalities  are  satisfied  for  all  s  in  [0, 1]: 


<*u(s,s0(t),<T(t),LT)  <  ol7(s,50(t),o(t),Z,r), 


(53) 


and 

Pu{s,si(.t),<r{t),LT)  <  Pu{s,si{t),cr{t),LT).  (54) 

Proof:  In  order  to  prove  (53)  and  (54)  under  the  conditions  of  Corollary  1,  it  suffices  to 
show  that  the  following  inequality  involving  the  Chernoff  distances  is  true 


Cj(a,ix)  >  Cj{s,Lt) 


(55) 


for  j  —  0, 1,  where  <?_,(«,  Lt)  are  the  Chernoff  distances  for  the  matched  case.  Indeed,  this 
is  satisfied  for  j  =  0, 1,  because  as  shown  in  the  proofs  of  Propositions  4  and  2: 

[T  -  ^o(0->2r/j\-|2  J,  ^  fT  r^l(0  ~  ®°(0  ■>  2r  -  (i \i2 

Jo {  >(oF"}  H  )]  dt  -  Jo  {~~[ [ 7(1)1 
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and 


«i(0  -  «o(<) 

l*(0)2  “ 


-  «o(0 

m? 


where  $o(0  =  so(t). 

Similarly  we  can  prove  the  above  proposition  under  the  conditions  of  Corollaries  2,  3, 
and  4.  The  following  proposition  holds  for  any  of  the  robust  tests  described  in  Corollaries 


1-4. 

Proposition  7:  For  problem  V\  with  the  uncertainty  models  for  the  signal  and  the  noise 
described  in  Corollaries  1-4,  if  the  likelihood  ratio  tests  of  (37)  or  (41)  are  employed,  then 
the  Chernoff  upper  bounds  on  the  error  probabilities  of  approach  zero  exponentially  with 
increasing  T  despite  the  uncertainty  . 

Proof:  Because  of  the  definitions  (45)-(46)  and  (55)  (which  was  essential  for  proving 
Proposition  6),  it  suffices  to  show  that 


Cj(s,  Lt)  >  0 


(56) 


for  all  s  in  [0, 1].  Indeed,  from  (50)  with  cl  and  C2  as  defined  there  but  with  Sj(t)  replaced 
by  Sj  for  j  =  0, 1  we  obtain 

Ca(s,tT)  =  1[ciJ2  +  c2«-sj f  3o(l)*] 


mop 

In  .  rT  [«.(*)-*>(*)]  V  n 
-  ?  ^  Vo  WW~ 


and  similarly  from  (51) 


i,  ,  rT ii(Q - s0(t) . 


Ci(s,Lt)  =  ^[cis2  -  c2s  +  s  J 


1  /,  x  fT  [h(t)  -  S0{t)]2 


5i(<)dt] 


/o  2  [d(t)]2 

Notice  that  the  Chernoff  distances  in  the  matched  case  are  strictly  nonzero  for  lj(t)  / 
s0(t),  0  <  t  <T.  Also,  C\(s,Lt)  =  Cq{s,Lt)  for  all  s  in  [0,1]. 
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If  the  limit  limr-^oo  j  Iq  exists,  then  Cj(a,L)  (j  =  0, 1)  becomes  inde¬ 

pendent  of  T  and  thus  the  Chernoff  bounds  on  the  error  probabilities  decrease  exponen¬ 
tially  to  zero  with  constant  rate,  as  T  increases. 
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IV.  MINIMAX  ROBUST  BLOCK  DISCRIMINATION  FOR  V2 

The  criteria  employed  for  the  robust  design  in  this  section  are  the  Chernoff  bounds 
on  the  error  probability  and  the  divergences  induced  by  them.  The  divergences  under 
mismatch  are  defined  as  follows: 

I0(L)  =  C'0(s,  £)|,=0  =  -  /  In  ^ dP0 ,  (57) 

J  dro 

and 

hit)  =  C[(s,L)\,= o  =  /  In  ^r-dPu  (58) 

J  dPo 

where  '  denotes  derivative  with  respect  to  a;  similarly  the  regular  (matched)  diver¬ 
gences  are  defined  as 

i0{L)  =  C'0{s ,  L) |,=0  =  -  /  In  ^dPo  >  0,  (59) 

J  dro 

and 

h(L)  =  C'(a,  i)|,=0  =  /in  dP1  >  0,  (60) 

J  dro 

The  strict  positivity  of  the  matched  divergences  for  P\  ^  Po  can  be  easily  established 
with  the  help  of  Jensen’s  inequality.  In  [10]  the  divergences  were  considered  in  the  context 
of  robust  discrimination  problems  from  discrete-time  observations.  The  usefulness  of  the 
divergence  distance  measures  lies  in  (i)  their  relation  to  the  Chernoff  bounds  on  the  error 
probabilities  of  the  hypothesis  testing  problems,  which  enables  us  to  carry  over  inequal¬ 
ities  valid  for  the  divergencxes  to  similar  inequalities  on  the  Chernoff  distances,  and  (ii) 
the  simplicity  of  their  form  which  enables  us  to  prove  dominance  properties  (inequalities 
for  minimax  robustness)  for  them,  which  are  impossible  to  prove  directly  for  the  error 
probabilities  or  the  Chernoff  bounds  on  them  for  a  variety  of  problems  with  statistical 
uncertainty. 
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From  the  form  of  the  Chernoff  bounds  (45)-(48),  we  notice  that  the  the  exponents  of 
the  bounds  are  convex  functions  of  s  in  [0, 1].  Therefore,  in  order  to  prove  inequalities  of 
the  form  (53)-(54)  which  are  equivalent  to  inequalities  of  the  form  (55) 

Cj(s,  £t)  >  £j(*»  It) 

Vs  €  [0,s]  for  some  s  >  0,  we  need  to  prove  the  equivalent  inequalities  between  the 
mismatch  and  matched  divergences: 

Ij(Lr)  =  Cj(s,  Lt)\,=o  >  6'(s,  Lr)Uo  =  Ij(LTy,  j  =  0, 1,  (61) 

Thus  the  final  criterion  turns  out  to  be  the  divergences  under  two  hypotheses,  which 
can  be  easily  obtained  from  the  likelihood  ratio  function  (14)  and  the  following  form  of 
solution  of  Yt  in  P2  (see  [4.4  of  9])  under  Hj  ( j  =  0,1) 

Yt  =  exp[  f  Si(t)dt]X o  +  /  exp[  /  «»(ri)dri]dWT,  0  <  t  <  T  (62) 

Jo  Jo  Jr 

Xq  has  a  standard  normal  distribution  and  Yt  is  independent  of  the  noise  NT  =  W (r)  for 
t  >t  under  either  hypothesis.  Substituting  (62)  into  (14)  and  using  the  independence  of 
Yt  and  WT  for  0  <  r  <  t,  we  have  the  following  general  forms  for  the  mismatch  divergences 
of  V-x  under  the  assumption  Ej\Y?)  <  oo  for  j  =  0, 1: 

Io(so(t),si{t),cr{t),LT)  =  -Eoilnix} 

=  -Eo<-  £  +  io(<)  -  2 Mt)]Y?4t 

+  rTm-mYtWiit) 

Jo  <r(t) 

=  £  my,2)* 
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(exp[2  JQ  MTi)dTi]  +  Jo  exp[2^  50(rj)<fri][ff(r2)]2<ir2}d< 

(63) 


/i(«o(0.«i(0»?(0  >Lt) 


Ei  {In  Lj } 

E' [  Jo  -  h(t)  -  km;,dt 

^  s:  m_^YtWtdt) 

JJ  ~ -  »■(')  - 

Jo  2[6(I)P  ^2>1^  "  J,W  -  *>W1 

{exp[2  jf  si(rl)dTi]  +  exp [2J  Si{Ti)dTi)[a{T2)fdT2}dt 

(64) 


where  .EjfU^VF,,]  =  ct(u)<t(v)6(u  —  v)  (j  =  0,1).  Furthermore,  the  matched  divergences 
are 


Io{h(t)>h{t),p(t),LT)  =  -EoQvLt} 

=  fT  Ml  ~  *>(<)]2 

/o  2[a(t)]2 

{exp[2  jf  s0(ri)dT1]  +  exp[2 ,J  5o(Ti)dr1][CT(T2)]2dr2}dt 

(65) 


Ji(^o(<))^i(0>&W>£t)  =  £i{lnZ:r} 

_  fT  (ji(0  -  *q(012 

A  2[<x(0]2 

{exp [2  jf  5i(r1)dr1]  + jf  exp[2^  Si(T1)dr1][d(r2)]2dr2}df 

(66) 
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Within  this  framework,  we  pursue  the  least-favorable  pair  within  the  class  of  likeli¬ 
hood  ratio  tests  for  the  divergences,  and  thus  for  the  Chernoff  bounds  on  the  two  error 
probabilities,  for  different  uncertainty  classes.  It  turns  out  that  we  can  prove  the  following 
proposition,  in  which  the  signal  So(t)  under  hypothesis  Hq  is  a  known  time  function  and 
the  noise  quantity  tr(<)  in  V2  belongs  to  the  class: 


El  =  MOKO  >  °l( <)»  €  [o,r]> 


(67) 


In  comparing  this  class  to  that  of  (11)  we  notice  that  it  is  characterized  by  a  minimum 
rather  than  a  maximum  element. 

Proposition  8:  For  problem  V2  with  signal  uncertainty  within  the  class  (10)  and  noise 
uncertainty  within  the  class  (67),  if  Jj(t)  =  Siti(t)  >  sQ(t)  =  s0(t)  and  a(t )  =  ox(t),  then 
the  likelihood  ratio  test  based  on 

= ■"  sr f  - f  (68) 

is  minmax  robust  for  the  divergences  and  thus  for  the  Chernoff  bounds  on  the  error 
probabilities  under  two  hypotheses;  that  is 


Ij(s0(t),  si(t)<r(t),  Lt)  >  si (<)»  ^(0>  ir);  j  =  0,1  (69) 

Proof:  By  comparing  (63)  with  (65)  and  (64)  with  (66),  respectively,  we  can  readily 
establish  the  desired  inequalities  (69)  for  j  =  0, 1. 

The  following  proposition  follows  directly  from  Proposition  8  in  the  same  way  that 
Proposition  7  followed  from  Proposition  6: 

Proposition  9:  Under  the  conditions  of  Proposition  8  the  Chernoff  upper  bounds  on  the 
error  probabilities  approach  zero  exponentially  with  increasing  T  despite  the  uncertainty. 
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V.  MINIMAX  ROBUST  SEQUENTIAL  DISCRIMINATION  FOR  HOMO¬ 
GENEOUS  Vi 

Here  we  design  a  minimax  robust  sequential  discriminator  based  on  a  robust  version 
of  the  sequential  probability  ratio  test  (SPRT)  for  the  case  that  the  signals  and 
the  variance  of  noise  in  V\  are  time  invariant,  that  is,  when  So(^)  =  ■®o>  si(t)  =  «i  and 
<r(t)  =  a\  in  this  case  the  observations  follow  the  equation 

dYt  =  Sidt  +  dWu  i  =  0, 1 


where  Wt  is  a  Wiener  process  with  zero  mean  and  intensity  function  a.  As  in  the  previous 
sections,  symbols  with  a  hat  denote  least-favorable  pairs  of  probabilities  and  the  likelihood 
ratios  based  on  them,  symbols  without  a  hat  denote  the  actual  operating  conditions  of 
V\.  From  Lemma  2  [eq.  (13)]  we  know  that  for  this  particular  case,  the  likelihood  ratio 
function  has  the  following  form 


In  Lt  =  In 


dPi 


dP0 

Jj  -  So 


<72 


fT  Ytdt  - 
Jo 


S2  -  s2 
si  so . 

2  a2 


(70) 


for  observations  over  the  interval  [0,  r].  In  this  section  we  also  assume,  without  loss  of 
generality,  that  Yo  —  0.  Let  a  and  /3  be  the  error  probabilities  under  hypotheses  Hq  and 
H i,  respectively,  and  A,B  be  the  two  thresholds  involved  in  the  SPRT  when  the  actual 
operating  conditions  are  matched  to  the  assumed  signals  and  noise  variance  Si,So)  and  o. 

The  SPRT  operates  as  follows:  if  Lr  >  B  hypothesis  H\  is  decided;  if  LT  <  A  hypoth¬ 
esis  Ho  is  favored;  finally,  if  A  <  LT  <  B  more  observations  are  collected.  The  time  r  that 
tany  of  the  two  thresholds  is  crossed  is  termed  the  stopping  time.  According  to  [11]  for 
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continuous-time  SPRTs,  a,  (3,  A  and  B  are  related  as  follows 


and 


A  = 


1  —  a 


As  usual,  we  make  the  regularity  condition  that 


a  -f  (i  <  1 


(71) 


(72) 


(73) 


i.e.,  the  false-alarm  probability  (a)  is  smaller  than  the  power  (1  —  /?)  of  the  test.  From 
(71)-(73)  one  obtains  that  0  <  A  <  1  <  B. 

Under  mismatch,  we  have  the  following  definitions  for  the  error  probabilities 

a  =  Po(LT>B),0  =  P1(Lr<A)} 

while  for  the  matched  case 

&  =  Po(Lr  >B)J=  Pl(Lr  <  A) 

Under  mismatch  the  expected  stopping  times  under  the  two  hypotheses  are  denoted 
by  Ej{r\LT }  ( j  =  0,1);  for  the  matched  case  they  are  denoted  by  Ej{r\LT}  ( j  =  0,1). 
Next  we  derive  the  following  proposition  for  the  situation  characterized  by  mismatch: 
Proposition  10:  Suppose  A  and  B  (0  <  A  <  1  <  B)  are  the  thresholds  of  the  SPRT  for 
a  homogeneous  V\.  If  the  signals  and  the  noise  variance  (ao,si,d)  are  used  in  determining 
the  likelihood  ratio  of  (70)  which  is  employed  by  the  SPRT,  while  the  actual  operating 
conditions  are  characterized  by  (so>si,o),  then  the  following  identities  hold 


E0{r\ir) 


—2  a7u(a,0,  a) 

(h  -  Jo)(2s0-ii  -  so)’ 


(74) 
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and 


where 


and 


pWf  2 0*0)0, a, &) 

Eli  \Er}-  ^-so^-h-ioY 

(75) 

_  f  x  _  A-H-l 

v  ’  T>  ~  (B~2A)^  -  1 

(76) 

/?(si,<7,  LT)  =  —  -  — . - - - 

v  '  (A-'B)*- 1 

(77) 

1  ~  X  X 

u>(x,y,x)  =  (1  -  i)ln — —  +  3  In- - r, 

y  i  -  y 

(78) 

cr2(2s0  -  Si  -  s0) 

70  “  -  so)  ’ 

(79) 

<72(2si  -  Si  -  So) 

7i  = - o T- — 

o\st  -  s0) 

(80) 

Proof:  First  we  prove  (74)  and  (75).  Let 

{1  if  t  <  t 
0  if  t  >  t 

be  the  indicator  function  of  the  set  {t  <  r},  then 


E0 


{IoYtdt]  = 


oij™  Ytltdt) 

=  r  E0{Yt}E0{It}dt 
Jo 

=  S0  f  Po{t  >  t)dt 
Jo 

=  SqE0{t } 


where  we  used  the  fact  that  Yt  and  It  are  independent  (this  is  a  result  of  the  independence 
of  Yt  and  Yu  for  (  /  «,  the  stopping  time  r  depends  on  all  Yu  for  0  <  «  <  r  but  not  on 
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the  current  value  Yt).  Then,  under  the  actual  operation  conditions,  we  have 


and 


(sx  -  So)(2s0  -  Si  -  So) 
2  a2 


Eq{t\Lt j 


JB0{ln~}  =  Po{Lt  >  B)\nB  +  P0{LT  <  A)\nA 
dPo 

=  a  In  *  .  -  4-  (1  -  a)  In 

a  1  —  a 

=  -w(a,4,o) 


Thus  (74)  follows  from  the  above  two  expressions  for  Eq{\rLt}.  Eq.  (75)  can  be  shown 
in  a  similar  way;  the  proof  is  omitted. 

Next  we  prove  (76)  by  using  the  method  provided  in  Section  6.4  of  [12]  and  ornit  the 
proof  of  (77).  Let  us  define  h\  =  (s*  -  s0)/a2  and  h2  =  (sj  -  si)/ 2cr2.  Then  under  the 
assumption  Yq  =  0,  we  have 


a(s0,<r;so, Si, a)  =  Po({LT  crosses  \nB  before  ln/4}|Yo  =  0) 

=  P0({Zt  crosses  ln5  before  \nA)\Zo:=0) 

where  Zt  =  hi  f£  Yudu  -  h2t  and,  under  Ho,  Zt  has  normal  distribution  with  mean  (/ii-so  — 
h2)t  and  variance  h\o2t.  Then  according  to  Section  6.4  of  [11]  we  can  think  of  a  as  a(0), 
where  o (z)  is  a  function  of  the  initial  condition  Z0  =  z  and  write  the  following  expression 


a(z)  =  E0[a(z  +  Z)]  +  o(h) 


where  o(h)  denotes  the  probability  that  the  process  would  have  already  crossed  one  of  the 
barriers  by  the  time  h  and  Z  =  Zh  —  Zq  =  Zh  -  z.  Notice  that  Z  has  a  normal  distribution 
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with  mean  (h\So  —  li2)h  —  z  and  variance  h\a2h.  Expanding  a(z  +  Z )  around  z  +  Z  —  z 
inside  the  expectation  Eq{}  in  the  previous  equation  we  obtain 

a(z)  =  Eo[a(z)  +  a'{z)Z  +  a"(z)Z2/ 2  +  •••]  +  o(h) 


where  a'(z),  a"(z)  •  ■  ■  denote  the  nth-derivative  (n  =  1,2,---)  of  a(z)  with  respect  to  z. 
Letting  h  —>0  and  taking  the  expectation,  we  have  the  following  differential  equation 


+  5 =  o 

2crz[Si  -  s0)  2 

By  solving  this  differential  equation  with  respect  to  a(z)  and  using  the  appropriate  bound¬ 
ary  conditions  for  the  two  barriers,  a;(lni?)  =  1,  a(ln  A)  =  0,  we  obtain  for  z  =  Zq  =  0 

,  x  -  1 

a  =  o(0)  =  - - - - 

v  J  (B-xAyi o  - 1 

with  70  defined  in  (79).  This  completes  the  proof  of  this  proposition. 

Equipped  with  the  above  proposition,  we  now  state  and  prove  the  main  result  of  this 
section  for  the  minimax  robust  discriminator  of  homogeneous  V\ . 

Proposition  11:  For  problem  V\  with  signal  and  noise  uncertainty  within  the  homoge¬ 
neous  classes  (9)-(  1 1 ) ,  suppose  si  =  s\l  >  sou  —  and  o  =  cq/,  and  that  the  (73)  and 
the  following  conditions  on  the  desired  error  probabilities  of  the  SPRT  are  satisfied: 


ln[(l  -  0)/&]  >  4’ln[(l  -  a)(l  -  $)/(&$))  (81) 

ln[(l  -  &)/$]  >  a  ln[(l  -  a)(l  -  $)/{&$)]  (82) 

then  the  SPRT  based  on  the  likelihood  ratio  of  (70)  is  minimax  robust  with  respect  to 

both  the  error  probabilities  and  the  expected  stopping  times  under  the  two  hypotheses; 

that  is 

1  —  A 

a(s0,(T,LT)  <  a(s0,cr,LT)  =  a  =  ~r - j,  (83) 
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(84) 


P(sUa,Lr)  <  P(suO,lr)  =  p  = 


Mb  - 1) 

B-A  ’ 


and 


Ej{T\lr)  <  i  =  0,l 


(85) 


Remark  1:  Conditions  (81)-(82)  are  not  too  restrictive,  they  are  easily  satisifed  if  both 
a  and  /3  are  smaller  than  10-2. 

Remark  2:  The  inequality  in  (85)  is  not  an  inequality  in  the  strict  sense;  it  is  however  a 
valid  inequality  for  all  practical  purposes  when  (81)  and  (82)  are  satisified. 

Remark  3:  As  d  — *•  0  and  (3  — » *  0  the  stopping  time  t  -*  oo  and  (85)  become  inequalities 
in  the  strict  sense: 


E0{r\LT}  = 


2 d2  In  A 


-2 d2  In  A 


(sj  -  S0)(2s0  -  Si  -  So)  (Si  -  So)2 


=  Eq{t\Lt} 


(86) 


_  r  lf  ,  2d2lnB  2d2 In B  ^  ,  /Q7, 

E\{t\Lt}  =  - - - rr  <  jz - =  Ex{T\Lr)  (87) 

(sj  -  s0)(2si  -  Si  -  s0)  («i  -  ■So)2 

The  above  expected  stopping  times  under  the  conditions  a  — *  0  and  $  — ♦  0  are  termed 
the  asymptotic  speeds  of  the  SPRT. 

Proof:  We  only  show  (83)  and  (85)  for  j  =  1;  a  similar  proof  can  be  obtained  for  (84) 
and  (85)  for  j  =  0.  First  we  prove  (83).  From  Proposition  10  we  notice  that  a  depends 
on  (so,<r)  only  through  70.  Taking  the  derivative  with  respect  to  70,  we  have 

- [(l_^>)lnA  +  (i^  -  l)lnJ5]. 

d7o  {(R-M)to-  l}2^  ' 

Let  G(jo)  =  (1  -  fC°)ln  A  +  (A10  -  l)lnR.  Then  the  sign  of  the  above  derivative  is 
determined  by  the  sign  of  <j(7o).  Since  the  derivative  of  G( 70)  with  respect  to  70  is 
nonpositive 

.MfpA  -  in  A  In  BiA™  -  &*)  <  0 
■  djo 
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=  In iln B(A'W  -  2^>)  <  0 


for  70  <  0  (In  A  <  0  <  In  B),  then  G(7o)  is  a  decreasing  function  of  70  for  70  <  0,  and 
thus  G(7o)  >  G(0)  =  0  for  70  <  0.  Consequently,  0(70)  is  an  increasing  function  of  70  for 
7o  <  0-  From  the  uncertainty  classes  (9)  and  (11)  in  the  homogeneous  case:  So  <  $ou  =  «o 
and  a  <  cry  =  d,  and  thus  we  can  easily  show  that 

7o<— 1 

where  the  equality  holds  when  (so,<r)  =  (so,d).  Therefore, 

a  =  a( 70)  <  or(-l)  =  a 

which  proves  (83).  Similarly,  to  prove  (84)  we  first  show  that  /?  is  a  decreasing  function 
of  71  and  that  71  >  1. 

Next  we  show  (85)  for  j  =  1.  Since  from  (84)  P  <  p  we  use  the  assumption  (81)  to 
obtain 


M(1  -  /?)/«]  >  /? ln[(l  -  6)(1  -  P)/(aP)\  >  /?ln[(l  -  d)(l  -  P)I{&P)] 


Then  by  using  this  inequality  and  the  form  of  Ei[r\LT]  from  Proposition  10,  we  have 

Ei[t\Lt]  = 


r. - -yP—- - + 

(si  -  so)(2si  -  si  -  so)  a  1  -  a 


2  d 


■In 


1  -P 


(si  -  so)(2sj  -  si  -  so)  & 

Using  (81)  directly  and  the  form  of  Ei[t\Lt ]  from  Proposition  10  for  the  matched  case, 
we  obtain 

~  m’-5 


(h  -  s0)2  & 
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Both  the  above  equations  are  approximations.  However,  as  noted  in  Remark  1  above, 
conditions  (81)-(82)  are  easy  to  satisfy,  and  this  guarantees  that  the  deviations  of  the 
numerators  in  both  of  the  above  equations  from  the  common  value  In  are  negligible. 
Since  the  denominators  are  expected  to  be  substantially  different  for  almost  all  elements 
«o  and  Sj  in  the  uncertainty  classes,  for  all  practical  purposes  the  above  approximations 
maintain  the  distinct  features  of  the  two  quantities  and  allow  us  to  compare  them  with  each 
other  without  any  measurable  loss  in  accuracy.  Indeed,  by  using  the  fact  that  si  >  sn,  =  Sj 
we  obtain 

E1[t\Lt]  <  Ei[t\Lt] 

which  completes  the  proof. 
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VI.  CONCLUSIONS 


In  this  paper  motivated  by  practical  situations  in  target  discrimination  we  consid¬ 
ered  binary  hypothesis  testing  problems  in  continuous-time  characterized  by  observations 
which  (i)  consist  of  distinct  signals  in  additive  white  Gaussian  noise  or  (ii)  are  the  output 
of  stochastic  dynamical  systems  driven  by  white  Gaussian  noise.  The  statistics  of  the 
observations  in  both  problems  were  only  partially  known.  The  goal  of  the  paper  was  to 
robustify  the  continuous-time  discrimination  tests  against  statistical  uncertainty  in  the 
observations.  In  particular,  the  signals  in  the  first  model,  the  parameters  of  the  dynam¬ 
ical  systems  in  the  second  model,  and  the  autocorrelation  functions  of  the  noise  in  both 
models  were  modeled  to  belong  to  one  of  the  following  distinct  uncertainty  classes:  (1) 
classes  determined  by  2-alternating  capacities  and  (2)  classes  with  minimum  or  maximum 
elements.  These  uncertainty  classes  include  many  popular  models  as  subcases  and  can 
model  effectively  many  practical  situations.  In  the  course  of  reviewing  the  theory  of  2- 
alternating  capacities  we  also  extended  Huber’s  results  on  minmax  robustness  within  such 
classes  to  more  general  functionals  than  the  average  risk. 

We  then  derived  minimax  robust  discrimination  tests  (a)  with  a  fixed  observation 
interval  and  (b)  sequential  tests.  The  likelihood  ratios  of  all  the  robust  tests  depended  on 
the  least-favorable  pairs  of  parameters  in  the  aforementioned  uncertainty  classes  and  were 
shown  to  have  an  acceptable  level  of  performance  despite  the  uncertainty.  For  the  robust 
tests  with  a  fixed  observation  interval  the  performance  measures  considered  for  problem 
(i)  were  the  actual  error  probabilities  under  the  two  hypotheses  and  the  Chernoff  upper 
bounds  on  them;  the  latter  were  shown  to  preserve  their  desirable  asymptotic  properties 
(in  the  sense  that  they  approach  zero  exponentially  as  the  observation  interval  increases) 
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>t  r  .  (  f 


in  the  presence  of  the  uncertainties.  For  problem  (ii)  we  provided  results  only  for  the 
Chemoff  bounds  and  the  associated  divergences.  Finally,  we  derived  the  robust  sequential 
test  for  problem  (i)  with  time-invariant  parameters;  in  this  case  we  used  as  performance 
measures  the  error  probabilities  and  the  average  required  length  of  the  observation  interval 
under  each  hypothesis. 
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